SARS-CoV-2 infection dynamics and genomic surveillance to detect variants in wastewater - a longitudinal study in Bengaluru, India.

Background: Environmental surveillance (ES) of a pathogen is crucial for understanding the community load of disease. As an early warning system, ES for SARS-CoV-2 has complemented routine diagnostic surveillance by capturing near real-time virus circulation at a population level. Methods: In this longitudinal study conducted between January 2022 and June 2022 in 28 sewershed sites in Bengaluru city (∼11 million inhabitants), we quantified weekly SARS-CoV-2 RNA concentrations to track infection dynamics and provide evidence of change in the relative abundance of emerging variants. Findings: We describe an early warning system using the exponentially weighted moving average control chart and demonstrate how SARS-CoV-2 RNA concentrations in wastewater correlated with clinically diagnosed new COVID-19 cases, with the trends appearing 8-14 days earlier in wastewater than in clinical data. This was further corroborated by showing that the estimated number of infections is strongly correlated with SARS-CoV-2 RNA copies detected in the wastewater. Using a deconvolution matrix, we detected emerging variants of concern up to two months earlier in wastewater samples. In addition, we found a huge diversity in variants detected in wastewater compared to clinical samples. The findings from this study have been discussed regularly with local authorities to inform policy-making decisions. Interpretation: Our study highlights that quantifying viral titre, correlating it with a known number of cases in the area, and combined with genomic surveillance helps in tracking variants of concern (VOC) over time and space, enabling timely and making informed policy decisions. Funding: This work has been supported by funding from the Rockefeller Foundation grant to National Centre for Biological Sciences, TIFR) and the Indian Council of Medical Research grant to (FI) Tata Institute for Genetics and Society and Tata Trusts.

Integrative genome-wide analysis reveals EIF3A as a key downstream regulator of translational repressor protein Musashi 2 (MSI2).

Musashi 2 (MSI2) is an RNA binding protein (RBP) that regulates asymmetric cell division and cell fate decisions in normal and cancer stem cells. MSI2 appears to repress translation by binding to 3' untranslated regions (3'UTRs) of mRNA, but the identity of functional targets remains unknown. Here, we used individual nucleotide resolution cross-linking and immunoprecipitation (iCLIP) to identify direct RNA binding partners of MSI2 and integrated these data with polysome profiling to obtain insights into MSI2 function. iCLIP revealed specific MSI2 binding to thousands of mRNAs largely in 3'UTRs, but translational differences were restricted to a small fraction of these transcripts, indicating that MSI2 regulation is not triggered by simple binding. Instead, the functional targets identified here were bound at higher density and contain more 'UAG' motifs compared to targets bound nonproductively. To further distinguish direct and indirect targets, MSI2 was acutely depleted. Surprisingly, only 50 transcripts were found to undergo translational induction on acute loss. Using complementary approaches, we determined eukaryotic translation initiation factor 3A (EIF3A) to be an immediate, direct target. We propose that MSI2 downregulation of EIF3A amplifies these effects on translation. Our results also underscore the challenges in defining functional targets of RBPs since mere binding does not imply a discernible functional interaction.

Cytokeratins 5 and 17 Maintain an Aggressive Epithelial State in Basal-Like Breast Cancer.

Basal-like breast cancers (BLBC) are the most common triple-negative subtype (hormone receptor and HER2 negative) with poor short-term disease outcome and are commonly identified by expression of basal cytokeratins (CK) 5 and 17. The goal of this study was to investigate whether CK5 and CK17 play a role in adverse behavior of BLBC cells. BLBC cell lines contain heterogeneous populations of cells expressing CK5, CK17, and the mesenchymal filament protein vimentin. Stable shRNA knockdown of either CK5 or CK17 compared with non-targeting control in BLBC cells was sufficient to promote an epithelial-mesenchymal transition (EMT) gene signature with loss of E-cadherin and an increase in vimentin expression. Relative to control cells, CK5 and CK17 knockdown cells acquired a more spindle-like morphology with increased cell scattering and were more invasive in vitro. However, CK5 or CK17 knockdown compared with control cells generated decreased lymph node and lung metastases in vivo. Loss of CK5 or CK17 moderately reduced the IC50 dose of doxorubicin in vitro and led to increased doxorubicin efficacy in vivo. Single-cell RNA-sequencing of BLBC patient-derived xenografts identified heterogeneous populations of CK5/CK17, vimentin, and dual basal CK/vimentin-positive cells that fell on an EMT spectrum of epithelial, mesenchymal, and intermediate, respectively, whereas knockdown of CK5 transitioned cells toward a more mesenchymal score. IMPLICATIONS: This study supports that basal CKs 5 and 17 contribute to the adverse behavior of BLBC cells and could be an untapped source of therapeutic vulnerability for this aggressive disease.

Estrogens and Progestins Cooperatively Shift Breast Cancer Cell Metabolism.

Metabolic reprogramming remains largely understudied in relation to hormones in estrogen receptor (ER) and progesterone receptor (PR) positive breast cancer. In this study, we investigated how estrogens, progestins, or the combination, impact metabolism in three ER and PR positive breast cancer cell lines. We measured metabolites in the treated cells using ultra-performance liquid chromatography coupled with mass spectrometry (UPLC-MS). Top metabolic processes upregulated with each treatment involved glucose metabolism, including Warburg effect/glycolysis, gluconeogenesis, and the pentose phosphate pathway. RNA-sequencing and pathway analysis on two of the cell lines treated with the same hormones, found estrogens target oncogenes, such as MYC and PI3K/AKT/mTOR that control tumor metabolism, while progestins increased genes associated with fatty acid metabolism, and the estrogen/progestin combination additionally increased glycolysis. Phenotypic analysis of cell energy metabolism found that glycolysis was the primary hormonal target, particularly for the progestin and estrogen-progestin combination. Transmission electron microscopy found that, compared to vehicle, estrogens elongated mitochondria, which was reversed by co-treatment with progestins. Progestins promoted lipid storage both alone and in combination with estrogen. These findings highlight the shift in breast cancer cell metabolism to a more glycolytic and lipogenic phenotype in response to combination hormone treatment, which may contribute to a more metabolically adaptive state for cell survival.

The genome trilogy of Anopheles stephensi, an urban malaria vector, reveals structure of a locus associated with adaptation to environmental heterogeneity.

Anopheles stephensi is the most menacing malaria vector to watch for in newly urbanising parts of the world. Its fitness is reported to be a direct consequence of the vector adapting to laying eggs in over-head water tanks with street-side water puddles polluted by oil and sewage. Large frequent inversions in the genome of malaria vectors are implicated in adaptation. We report the genome assembly of a strain of An. stephensi of the type-form, collected from a construction site from Chennai (IndCh) in 2016. The genome reported here with a L50 of 4, completes the trilogy of high-resolution genomes of strains with respect to a 16.5 Mbp 2Rb genotype in An. stephensi known to be associated with adaptation to environmental heterogeneity. Unlike the reported genomes of two other strains, STE2 (2R+b/2Rb) and UCI (2Rb/2Rb), IndCh is found to be homozygous for the standard form (2R+b/2R+b). Comparative genome analysis revealed base-level details of the breakpoints and allowed extraction of 22,650 segregating SNPs for typing this inversion in populations. Whole genome sequencing of 82 individual mosquitoes from diverse geographical locations reveal that one third of both wild and laboratory populations maintain the heterozygous genotype of 2Rb. The large number of SNPs can be tailored to 1740 exonic SNPs enabling genotyping directly from transcriptome sequencing. The genome trilogy approach accelerated the study of fine structure and typing of an important inversion in An. stephensi, putting the genome resources for this understudied species on par with the extensively studied malaria vector, Anopheles gambiae. We argue that the IndCh genome is relevant for field translation work compared to those reported earlier by showing that individuals from diverse geographical locations cluster with IndCh, pointing to significant convergence resulting from travel and commerce between cities, perhaps, contributing to the survival of the fittest strain.

Neoadjuvant endocrine therapy expands stromal populations that predict poor prognosis in estrogen receptor-positive breast cancer.

The tumor microenvironment (TME) is an important modulator of response and resistance to endocrine therapy in estrogen receptor alpha (ER) positive breast cancer. Endocrine therapy is highly effective at reducing tumor burden and preventing recurrence in most estrogen receptor alpha (ER) positive breast cancers. Existing drugs work either directly by targeting tumor-cell ER or indirectly by inhibiting estrogen production in stromal cells with aromatase inhibitors (AI). However, many stromal cells also express ER and the direct impact of endocrine therapies on ER + stromal cells remain unclear. In this study, we investigated how neoadjuvant endocrine therapy (NET) directly effects stromal cells by measuring changes in stomal components of the TME that favor tumor progression. We previously defined two major subsets of tumor-associated stromal cells (TASCs): CD146 positive/CDCP1 negative (TASCCD146 ), CD146 negative/CDCP1 positive (TASCCDCP1 ), and generated a differentially expressed genes list associated with each type. Here, we applied the TASC gene list for classification and an algorithm that estimates immune cell abundance (TIMEx) to METABRIC transcriptomic data for ER + breast cancer patients coupled with multiplex imaging and analysis of paired tissue samples pre- and post- NET with the AI exemestane. TASCCDCP1 composition predicted for decreased patient survival in the METABRIC cohort. Exemestane treatment significantly increased expression of TASCCDCP1 and decreased expression of TASCCD146 . The posttreatment shift toward TASCCDCP1 composition correlated with increased macrophage infiltration and increased CD8+ T-cell, B cell, and general stromal components. The effectiveness of NET is currently based solely on the reduction of ER+ breast cancer cells. Here, we show NET displays clear TME effects that promote the expansion of the less favorable TASCCDCP1 population which are correlated with TME remodeling and reshaping immune infiltration supportive of tumor progression. Our findings highlight the need to further understand the role of endocrine therapy on TME remodeling, tumor progression, and patient outcomes.

Hidden genomic features of an invasive malaria vector, Anopheles stephensi, revealed by a chromosome-level genome assembly.

BACKGROUND: The mosquito Anopheles stephensi is a vector of urban malaria in Asia that recently invaded Africa. Studying the genetic basis of vectorial capacity and engineering genetic interventions are both impeded by limitations of a vector's genome assembly. The existing assemblies of An. stephensi are draft-quality and contain thousands of sequence gaps, potentially missing genetic elements important for its biology and evolution. RESULTS: To access previously intractable genomic regions, we generated a reference-grade genome assembly and full transcript annotations that achieve a new standard for reference genomes of disease vectors. Here, we report novel species-specific transposable element (TE) families and insertions in functional genetic elements, demonstrating the widespread role of TEs in genome evolution and phenotypic variation. We discovered 29 previously hidden members of insecticide resistance genes, uncovering new candidate genetic elements for the widespread insecticide resistance observed in An. stephensi. We identified 2.4 Mb of the Y chromosome and seven new male-linked gene candidates, representing the most extensive coverage of the Y chromosome in any mosquito. By tracking full-length mRNA for > 15 days following blood feeding, we discover distinct roles of previously uncharacterized genes in blood metabolism and female reproduction. The Y-linked heterochromatin landscape reveals extensive accumulation of long-terminal repeat retrotransposons throughout the evolution and degeneration of this chromosome. Finally, we identify a novel Y-linked putative transcription factor that is expressed constitutively throughout male development and adulthood, suggesting an important role. CONCLUSION: Collectively, these results and resources underscore the significance of previously hidden genomic elements in the biology of malaria mosquitoes and will accelerate the development of genetic control strategies of malaria transmission.

Insulin receptor substrate-1 (IRS-1) mediates progesterone receptor-driven stemness and endocrine resistance in oestrogen receptor+ breast cancer.

BACKGROUND: Progesterone receptors (PR) are potent modifiers of endocrine responses. In aberrant signalling cancer contexts, phosphorylation events dramatically alter steroid hormone receptor action. METHODS: The transcriptomes of primary tumours and metastases in mice harbouring ER+ breast cancer patient-derived xenografts (PDXs) were analysed following single-cell RNAseq. In vitro assays were employed to delineate mechanisms of endocrine resistance and stemness. RESULTS: A 16-gene phospho-Ser294 PR (p-PR) signature predicted poor outcome in ER+ breast cancer. Relative to primary PDX tumours, metastatic lesions expressed abundant p-PR and exhibited an activated PR gene programme with elevated expression of PGR and IRS-1. Breast cancer models of activated PR lost the expression of IGF1R and acquired insulin hypersensitivity with tamoxifen insensitivity. Activated p-PR+ breast cancer cells formed increased tumourspheres with enlarged ALDH+ and CD24-/CD44 populations. E2 induced PR/IRS-1 interaction and exchange of IGF1Rß for IRS-1 in p-PR-containing transcriptional complexes. Inhibition of IRS-1 or IR and inducible IRS-1 knockdown reduced tumourspheres. Endocrine-resistant models of luminal B breast cancer induced p-PR in 3D cultures and required PR and IRS-1 for tumoursphere formation. CONCLUSIONS: Phospho-PR-B cooperates with IRS-1 to promote outgrowth of endocrine-resistant and stem-like breast cancer cells. Targeting phospho-PR/IRS-1 crosstalk may block the emergence of endocrine resistance.

A Near-Chromosome Level Genome Assembly of Anopheles stephensi.

Malaria remains a major healthcare risk to growing economies like India, and a chromosome-level reference genome of Anopheles stephensi is critical for successful vector management and understanding of vector evolution using comparative genomics. We report chromosome-level assemblies of an Indian strain, STE2, and a Pakistani strain SDA-500 by combining draft genomes of the two strains using a homology-based iterative approach. The resulting assembly IndV3/PakV3 with L50 of 9/12 and N50 6.3/6.9 Mb had scaffolds long enough for building 90% of the euchromatic regions of the three chromosomes, IndV3s/PakV3s, using low-resolution physical markers and enabled the generation of the next version of genome assemblies, IndV4/PakV4, using HiC data. We have validated these assemblies using contact maps against publicly available HiC raw data from two strains including STE2 and another lab strain of An. stephensi from UCI and compare the quality of the assemblies with other assemblies made available as preprints since the submission of the manuscript. We show that the IndV3s and IndV4 assemblies are sensitive in identifying a homozygous 2Rb inversion in the UCI strain and a 2Rb polymorphism in the STE2 strain. Multiple tandem copies of CYP6a14, 4c1, and 4c21 genes, implicated in insecticide resistance, lie within this inversion locus. Comparison of assembled genomes suggests a variation of 1 in 81 positions between the UCI and STE2 lab strains, 1 in 82 between SDA-500 and UCI strain, and 1 in 113 between SDA-500 and STE2 strains of An. stephensi, which are closer than 1 in 68 variations among individuals from two other lab strains sequenced and reported here. Based on the developmental transcriptome and orthology of all the 54 olfactory receptors (ORs) to those of other Anopheles species, we identify an OR with the potential for host recognition in the genus Anopheles. A comparative analysis of An. stephensi genomes with the completed genomes of a few other Anopheles species suggests limited inter-chromosomal gene flow and loss of synteny within chromosomal arms even among the closely related species.

Analysis of circulating breast cancer cell heterogeneity and interactions with peripheral blood mononuclear cells.

For solid tumors, extravasation of cancer cells and their survival in circulation represents a critical stage of the metastatic process that lacks complete understanding. Gaining insight into interactions between circulating tumor cells (CTCs) and other peripheral blood mononuclear cells (PBMCs) may provide valuable prognostic information. The purpose of this study was to use single-cell RNA-sequencing (scRNA-seq) of liquid biopsies from breast cancer patients to begin defining intravascular interactions. We captured CTCs from the peripheral blood of breast cancer patients using size-exclusion membranes followed by scRNA-seq of enriched CTCs and carry-over PBMCs. Transcriptome analysis identified two populations of CTCs: one enriched for transcripts indicative of estrogen responsiveness and increased proliferation and another enriched for transcripts characteristic of reduced proliferation and epithelial-mesenchymal transition (EMT). We applied interactome and pathway analysis to determine interactions between CTCs and other captured cells. Our analysis predicted for enhanced immune evasion in the CTC population with EMT characteristics. In addition, PD-1/PD-L1 pathway activation and T cell exhaustion were predicted in T cells isolated from breast cancer patients compared with normal T cells. We conclude that scRNA-seq of breast cancer CTCs generally stratifies them into two types based on their proliferative and epithelial state and differential potential to interact with PBMCs. Better understanding of CTC subtypes and their intravascular interactions may help design treatments directed against CTCs with high metastatic and immune-evasive competence.

New generation breast cancer cell lines developed from patient-derived xenografts.

BACKGROUND: Breast cancer is a highly heterogeneous disease characterized by multiple histologic and molecular subtypes. While a myriad of breast cancer cell lines have been developed over the past 60 years, estrogen receptor alpha (ER)+ disease and some mutations associated with this subtype remain underrepresented. Here we describe six breast cancer cell lines derived from patient-derived xenografts (PDX) and their general characteristics. METHODS: Established breast cancer PDX were processed into cell suspensions and placed into standard 2D cell culture; six emerged into long-term passageable cell lines. Cell lines were assessed for protein expression of common luminal, basal, and mesenchymal markers, growth assessed in response to estrogens and endocrine therapies, and RNA-seq and oncogenomics testing performed to compare relative transcript levels and identify putative oncogenic drivers. RESULTS: Three cell lines express ER and two are also progesterone receptor (PR) positive; PAM50 subtyping identified one line as luminal A. One of the ER+PR+ lines harbors a D538G mutation in the gene for ER (ESR1), providing a natural model that contains this endocrine-resistant genotype. The third ER+PR-/low cell line has mucinous features, a rare histologic type of breast cancer. The three other lines are ER- and represent two basal-like and a mixed ductal/lobular breast cancer. The cell lines show varied responses to tamoxifen and fulvestrant, and three were demonstrated to regrow tumors in vivo. RNA sequencing confirms all cell lines are human and epithelial. Targeted oncogenomics testing confirmed the noted ESR1 mutation in addition to other mutations (i.e., PIK3CA, BRCA2, CCND1, NF1, TP53, MYC) and amplifications (i.e., FGFR1, FGFR3) frequently found in breast cancers. CONCLUSIONS: These new generation breast cancer cell lines add to the existing repository of breast cancer models, increase the number of ER+ lines, and provide a resource that can be genetically modified for studying several important clinical breast cancer features.

Cytokeratin 5 alters ß-catenin dynamics in breast cancer cells.

Estrogen receptor (ER) positive breast cancers often contain subpopulations of cells that express the intermediate filament protein cytokeratin 5 (CK5). CK5+ cells are enriched in cancer stem cell (CSC) properties, can be induced by progestins, and predict poor prognosis in ER+ breast cancer. We established through CK5 knockout and overexpression in ER+ breast cancer cell lines that CK5 is important for tumorsphere formation, prompting us to speculate that CK5 has regulatory activity in CSCs. To interrogate CK5 interacting proteins that may be functionally cooperative, we performed immunoprecipitation-mass spectrometry for CK5 in ER+ breast cancer cells. Focusing on proteins with signaling activity, we identified ß-catenin, a key transcription factor of the Wnt signaling pathway and cell adhesion molecule, as a CK5 interactor, which we confirmed by co-immunoprecipitation in several breast cancer models. We interrogated the dual functions of ß-catenin in relation to CK5. Knockout or knockdown of CK5 ablated ß-catenin transcriptional activity in response to progestins and Wnt stimuli. Conversely, CK5 induced by progestins or overexpression was sufficient to promote the loss of ß-catenin at the cell membrane and total E-cadherin loss. A breast cancer patient-derived xenograft showed similar loss of membrane ß-catenin and E-cadherin in CK5+ but not intratumoral CK5- cells and single-cell RNA sequencing found the top enriched pathways in the CK5+ cell cluster were cell junction remodeling and signaling. This report highlights that CK5 actively remodels cell morphology and that blockade of CK5-ß-catenin interaction may reverse the detrimental properties of CK5+ breast cancer cells.

Degenerate minigene library analysis enables identification of altered branch point utilization by mutant splicing factor 3B1 (SF3B1).

Cancer-associated mutations of the core splicing factor 3 B1 (SF3B1) result in selection of novel 3' splice sites (3'SS), but precise molecular mechanisms of oncogenesis remain unclear. SF3B1 stabilizes the interaction between U2 snRNP and branch point (BP) on the pre-mRNA. It has hence been speculated that a change in BP selection is the basis for novel 3'SS selection. Direct quantitative determination of BP utilization is however technically challenging. To define BP utilization by SF3B1-mutant spliceosomes, we used an overexpression approach in human cells as well as a complementary strategy using isogenic murine embryonic stem cells with monoallelic K700E mutations constructed via CRISPR/Cas9-based genome editing and a dual vector homology-directed repair methodology. A synthetic minigene library with degenerate regions in 3' intronic regions (3.4 million individual minigenes) was used to compare BP usage of SF3B1K700E and SF3B1WT. Using this model, we show that SF3B1K700E spliceosomes utilize non-canonical sequence variants (at position -1 relative to BP adenosine) more frequently than wild-type spliceosomes. These predictions were confirmed using minigene splicing assays. Our results suggest a model of BP utilization by mutant SF3B1 wherein it is able to utilize non-consensus alternative BP sequences by stabilizing weaker U2-BP interactions.

The Development and Use of Scalable Systems for Studying Aberrant Splicing in SF3B1-Mutant CLL.

Mutational landscape of CLL is now known to include recurrent non-synonymous mutations in SF3B1, a core splicing factor. About 5-10% of newly diagnosed CLL harbor these mutations which are typically limited to HEAT domains in the carboxyl-terminus of the protein. Importantly, the mutations are not specific to CLL but also present in several unrelated clonal disorders. Analysis of patient samples and cell lines has shown the primary splicing aberration in SF3B1-mutant cells to the use of novel or "cryptic" 3' splice sites (3SS). Advances in genome-editing and next-generation sequencing (NGS) have allowed development of isogenic models and detailed analysis of changes to the transcriptome with relative ease. In this manuscript, we focus on two relevant methods to study splicing factor mutations in CLL: development of isogenic scalable cell lines and informatics analysis of RNA-Seq datasets.